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Abstract 

LHC searches for non-standard Higgs bosons decaying into tan lepton pairs constitute a sensitive exper¬ 
imental probe for physics beyond the Standard Model (BSM), such as Supersymmetry (SUSY). Recently, 
the limits obtained from these searches have been presented by the CMS collaboration in a nearly model- 
independent fashion — as a narrow resonance model — based on the full 8 TeV dataset. In addition to 
publishing a 95% C.L. exclusion limit, the full likelihood information for the narrow resonance model has 
been released. This provides valuable information that can be incorporated into global BSM fits. We present 
a simple algorithm that maps an arbitrary model with multiple neutral Higgs bosons onto the narrow res¬ 
onance model and derives the corresponding value for the exclusion likelihood from the CMS search. This 
procedure has been implemented into the public computer code HiggsBounds (version 4.2.0 and higher). 
We validate our implementation by cross-checking against the official CMS exclusion contours in three 
Higgs benchmark scenarios in the Minimal Supersymmetric Standard Model (MSSM), and find very good 
agreement. Going beyond validation, we discuss the combined constraints of the rr search and the rate 
measurements of the SM-like Higgs at 125 GeV in a recently proposed MSSM benchmark scenario, where 
the lightest Higgs boson obtains SM-like couplings independently of the decoupling of the heavier Higgs 
states. Technical details for how to access the likelihood information within HiggsBounds are given in the 
appendix. The program is available at http://higgsbounds.hepforge.org 
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I. INTRODUCTION 

The search for Higgs bosons [T] continues to be a cornerstone of the physics program at the Large 
Hadron Collider (LHC). After the discovery of a Higgs boson by ATLAS [2] and CMS [3] it is crucial 
to find out whether the detected particle is part of a Higgs sector that contains several physical 
states. Higgs sectors of this kind are predicted in many theories of physics beyond the Standard 
Model (SM). For the understanding of the mechanism of electroweak symmetry breaking two 
complementary experimental endeavors are important: On the one hand the precise determination 
of the properties of the Higgs signal detected at around 125 GeV, and on the other hand the search 
for additional Higgs bosons. Both are crucial in the quest to identify the underlying physics. 
The existing limits from the Higgs searches at LEP, the Tevatron and the LHC already put very 
important constraints on the parameter spaces of different models that provide a Higgs-like state 
compatible with the detected signal. More data on both the detected signal and on searches for 
additional Higgs bosons will further enhance the sensitivity for discriminating possible scenarios of 
new physics from the SM and from each other. 

In order to facilitate the available experimental information from the Higgs searches at LEP, the 
Tevatron and the LHC, expressed in terms of relatively model-independent cross-section limits for 
testing a wide variety of theoretical models, the program HiggsBounds [IHZ! has been developed. 
Experimental information on the Higgs signal detected at a mass value of around 125 GeV is utilized 
in the sister program HiggsSignals [8] for testing the theoretical predictions from any kind of 
Higgs sector. The experimental information on the detected signal incorporated in HiggsSignals is 
turned into a likelihood, which is suitable for the inclusion into global fits (see, e.g., Refs. IMlj), 
where in addition many other observables are taken into account. In contrast, exclusion limits have 
traditionally been presented in terms of 95% C.L. limits, which a priori only provide the information 
whether a particular parameter point is excluded or not at the 95% C.L. by the considered search 
channel. In a global fit, where the predictions of a model are confronted with a large number of 
observables, it would usually be too restrictive to disregard a certain parameter point just because 
it falls outside of the 95% C.L. region of a single search channel. In fact, testing a large variety 
of observables one would expect that the measured values of some observables lie outside of the 
respective 95% C.L. regions for purely statistical reasons. It would therefore be very desirable if 
also negative experimental outcomes from Higgs searches were provided in terms of the likelihood 
information in the relevant parameters, instead of a simple binary rejection or acceptance at a 
certain confidence level (C.L.). Up to now, likelihood information was available in HiggsBounds 
only for the results from the LEP Higgs searches [7], while for all search channels at the Tevatron 
and LHC only 95% C.L. limits have been accessible. We report here on significant progress in this 
direction for the LHC Higgs boson search in the final state, which plays a central role in the 

search for additional Higgs bosons. 

Many models that can accommodate a SM-like Higgs boson at 125 GeV, such as the Mini¬ 
mal Supersymmetric Standard Model (MSSM) or the various types of Two-Higgs-Doublet Models 
(2HDM), predict additional Higgs bosons that decay predominantly into SM fermions. Therefore, 
LHG searches for new neutral Higgs bosons decaying to play a crucial role. In particular 
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within the MSSM, these searches lead to large excluded regions in the parameter space. The high¬ 
est experimental sensitivity occurs for smaller values of the CP-odd Higgs boson mass, Ma, and 
larger values of tan/3, the ratio of the two vacuum expectation values Hade]. 

One complication that arises for this search channel is the fact that two different production 
modes, gluon fusion and b quark associated production, can both be important. Their individual 
contributions to the signal rate can vary strongly over the parameter space. Since the acceptances 
of these two channels can also be very different, a two-dimensional cross section interpretation 
for the T~^T~ final state is desirable as a basis for (close to) model-independent exclusion limits. 
Recently, the CMS collaboration has published the likelihood information for their Higgs boson 
search in the final state m- The likelihood is given as a function of the two relevant Higgs 

production channels, gluon fusion and b quark associated production, for various mass values of 
the narrow resonance assumed for the signal model. 

In this paper we investigate the application of this new experimental information for testing the 
theoretical predictions of extended Higgs sectors and its incorporation in global fits. We develop a 
simple algorithm that maps an arbitrary model with in general several neutral Higgs bosons onto 
the narrow resonance model. In this way the corresponding value of the exclusion likelihood from 
the CMS search for the tested model can be determined. We furthermore describe the inclusion of 
this likelihood information into the publicly available Fortran code HiggsBounds [IHIj- For nearly 
any model under consideration, HiggsBounds provides an evaluation of the exclusion likelihood 
for a model parameter point based on the information from Ref. |15) . While the new likelihood 
information goes well beyond the standard test whether a particular parameter point is excluded 
at the 95% C.L., the likelihood information can also be employed to run HiggsBounds in this 
“standard” mode. In this case, HiggsBounds determines the parameter region that is excluded at 
the 95% C.L. based on all available searches, including the new t~^t~ result from CMS. The new 
version of HiggsBounds can be used together with its sister program HiggsSignals [8] in order 
to take into account both the information from search limits and from the detected signal for a 
comprehensive test of Higgs phenomenology. Both codes are available at: 


http://higgsbounds.hepf orge.org 

The paper is organized as follows. In Section [IT] we briefly summarize the experimental results 
that are used as input for our investigations. Details of the employed algorithm and the imple¬ 
mentation of the exclusion likelihood of Ref. [I5] into HiggsBounds are given in Section Hill The 
validation in various MSSM Higgs benchmark scenarios is discussed in Section HVl As an appli¬ 
cation, in Section m we investigate the constraints on a certain benchmark scenario in the MSSM 
that are obtained from using the new exclusion likelihood in combination with the information on 
the detected signal incorporated in HiggsSignals. We conclude in Section IVTl Finally, all relevant 
information needed to run HiggsBounds to obtain the likelihood information for the Higgs 

search channel for any parameter point under investigation are contained in an Appendix. 
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II. EXPERIMENTAL RESULTS 

This section briefly summarizes the experimental results from the CMS non-standard Higgs 
search in the rr final state m, that we have used as starting point for our investigation and that 
we have implemented in HiggsBounds. The search analysis of CMS is carried out in two separate 
selection categories: One requiring the presence of at least one 6-tagged jet, and one without the 
presence of a 6-tag. The former category is enriched by the production of a Higgs boson, denoted 
generically by (f), in association with two 6 quarks, gg —)• bb(j), while the latter is dominated by the 
gluon fusion process, gg —?■ (j). Hence, the search features sensitivity to the two different production 
modes separately, which enables the presentation of the search results in terms of individual signal 
strengths in both production modes for all tested Higgs boson masses. Separate information on 
the two production modes is an indispensable ingredient for enabling the presentation of (close to) 
model-independent exclusion limits or measurements, in case of a discovery. 

The data is further classified in categories defined by the two r lepton decay modes: er/j, gTh, 
eg, gg and ThTh, where Th denotes a hadronically decaying r lepton. Using a maximum likelihood 
technique, an estimator for the true rr invariant mass, rrirr, is reconstructed from the momenta 
of the visible r decay products and the missing transverse energy in the event. The uncertainty of 
the nirT reconstruction is estimated to be around 20% when averaged over all decay modes m- 

The resulting uirT spectrum in all categories (6-tag and r decay) separately is then subject to a 
profile likelihood analysis m, where the background parametrization, obtained from control region 
data and Monte Carlo simulation, and the signal shape parametrization are fitted simultaneously 
to the reconstructed mass spectrum. The fit is performed individually for test masses between 
90 GeV and 1 TeV, and the results are interpolated between the test masses. 

CMS interprets the results both in a nearly model-independent waj|^ for a single narrow res¬ 
onance (/>, and, in a model-specific context, in the MSSM, where three neutral Higgs bosons h, 
H and A potentially comprise the signal. The latter interpretation performs a likelihood ratio 
hypothesis test for the two hypotheses of a single SM-like Higgs boson at rrih = 125 GeV with 
exact SM properties and, alternatively, for the signal consisting of all three neutral Higgs bosons 
of the MSSM. In the latter case, the mrT distributions of the h/H/A —)■ rr decays are combined 
before the calculation of the likelihood. Note that, whereas the model-specific limits for the MSSM 
are based on the full integrated luminosity of the combined 7-1-8 TeV dataset, the results for the 
single narrow resonance model are obtained from only the 8 TeV dataset. 

Since HiggsBounds is designed to test any extended Higgs sector, with any coupling properties 
of the 125 GeV Higgs candidate (if the model under consideration provides such a candidate; 
obviously the phenomenological interest in other models is rather limited) and any masses and 
properties of the remaining Higgs spectrum, the nearly model-independent single resonance results 
are chosen for the implementation in HiggsBounds. On the one hand, this is the only possibility 
unless one is willing to adopt further model-dependent assumptions. On the other hand, it will in 

^ The presentation of the search results in terms of a limit on the inclusive total cross section times branching ratio 
inevitably involves a slight model dependence from the extrapolation to the inclusive quantity. In other words, 
the expectation of the kinematic distributions of the signal and/or background is model dependent. 
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general yield weaker, i.e. more conservative, limits than a dedicated analysis taking into account 
the full structure of the considered model. For example, in the MSSM this will yield a conservative 
limit whenever either mu/A ~ or, more generally, whenever the model predicts that more than 
one Higgs boson have a non-negligible signal yield and contribute in different regions in rrirr- Since 
the likelihood is constructed for single resonances, such a case cannot be properly reconstructed 
from the likelihood. However, if e.g. the heavy MSSM Higgs bosons H and A contribute at the 
same point in rurr^ their signal rates can be added before interpreting the likelihood. In this case 
the implementation is not necessarily conservative. A detailed study on the applicability of these 


limits to the MSSM benchmarks is presented in Section IV 


The profile likelihood analysis follows the standard implementation: The test statistics is given 


by 


^ C{N\(i- s{m) + b,e) 

Here N is the number of measured events, h and s{m) the number of expected background and 
signal events for a given resonance mass hypothesis m, the signal strength modifier, and 9 are 
the nuisance parameters decribing the systematic uncertainties. 9^ maximizes the likelihood in 
the numerator given a certain value of /r, whereas the likelihood reaches its global maximum at jl 
and 9, which is given in the denominator. The constraint 0 < fi < fj, is employed to not penalize 
the model for a possible excess of the data over the signal plus background prediction. It should 
be noted that the signal yield contains two independent components, corresponding to the two 
production modes gg ^ (j) and gg —)■ bbcj). Thus, g and s(m) are two-component vectors. 

Using toy Monte Carlo techniques or asymptotic expressions for large statistics, the expected 
probability distributions P(q^|hypothesis) can be constructed for the test statistics given above. 
In the model-independent analysis used here, the hypothesis either consists of Pi = g- s{m) -|-b for 
the case of the presence of a single narrow resonance with a given mass m and signal yield g ■ s{m), 
and of H 2 = b for the background. Using these hypothesis definitions and the observed value of 
the likelihood ratio, given by Eq. a with N given by the actual observed number of events, 
N = Aobs) the likelihood ratio technique can be used to define the CL^ as 


_ P{qn > s{m)+ b) 

P{qn > 


( 2 ) 


In a stand-alone search, the criterion CL^ < a is then used to exclude the presence of a signal at 
1 — a confidence level (C.L.). For the expected limit, the observed data in the calculation of 
is replaced by the median of the background-only expectation for q^. A model-independent limit, 
e.g. at the 95% C.L., can then be derived by varying g until CLg = 0.05. The value g at which 
this happens then represents the signal strength modiher which is just allowed at the 95% C.L.. 

For any given model, HiggsBounds reconstructs the predicted signal yield s(m) from the the¬ 
oretical input provided by the user, and obtains the corresponding value of the test statistics 
q^j_=i (or simply denoted ^modei) from the CMS likelihood data. The details of this procedure will 
be described in the following section. Note that HiggsBounds directly employs the expectation 
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and observation of the test statistics, and respectively, as the provided CMS data does 
not allow for a full reconstruction of the CL^ value. Nevertheless, in the limit of large numbers, 
the test statistics can be approximated by a chi-squared differences function above minimum, 
Q/i ~ = x^ - Xmin> as 

-2\nC{N\i^- s{m) + b,§^), (3) 

xLn ~ - 21n£(iV|//• s(m)-h (4) 

Thus, for example, we approximately obtain the two-dimensional limit at the 68% and 95% C.L. 
when the test statistics takes the values 2.28 and 5.99, respectively. 

As an example, we show the observed likelihood distribution, for test masses of 125 GeV 
and 300 GeV in Fig. [Tj We also indicate the approximated 68% G.L. and 95% G.L. limits by 
contour lines. It should be kept in mind that these contours are only for illustrational purposes. In 
HiggsBounds the full likelihood information is used, and the specihc limit at a certain C.L. can 
be easily obtained from this information. 

The implementation of the likelihood for the CMS rr search differs in two signihcant ways 
from the implementation of the LEP Higgs search x^ in HiggsBounds, which is already available 
since version 2.0.0 [5]. In the LEP implementation, each Higgs search channel, comprised of one 
production mode and one Higgs boson decay mode, is treated separately, thus no combination of 
production modes is applied or possible for the user. In addition, in the LEP implementation the 
X^ is estimated from the CL^+ft value in each channel at the given signal strength prediction using 
Gaussian approximations. In contrast, for the CMS </>—?• rr search the exact values of the test 
statistics as presented by CMS are used and properly combined for both production modes. 


III. LIKELIHOOD RECONSTRUCTION FOR EXTENDED HIGGS SECTORS 

Eor the construction of the exclusion likelihood from the iL —)• rr search, we make use of 
the following quantities: Eor each neutral Higgs boson, hi {i = l,...,iV), in a model with N 
neutral Higgs bosons, we have a prediction of the mass, m* (where the relevant range is currently 
mi G [90, 1000] GeV), the gluon fusion production cross section, a{gg —>■ hi), the cross section for 
production in association with b quarks, cr{gg —>■ bbhi), and the branching fraction — )• rr). 

The main algorithm for the likelihood reconstruction proceeds as follows: 

1. Signal rates of multiple Higgs bosons that cannot be resolved by the experimental analysis 
are added. We thus combine the signal predictions for two Higgs boson hi and hj (j / f), if 

\mi — mj\ < 20% • max{mi,mj). (5) 

Each Higgs boson can appear in different such combinations. Eor each combination, also 
called Higgs cluster and labeled with capital characters in the following, we evaluate the 
physical quantities as follows: We assume that the total rates are given by the incoherent 
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(a) Observed exclusion likelihood, and 
approximated exclusion contours for a resonance mass 
of = 125 GeV. 


(b) Observed exclusion likelihood, and 
approximated exclusion contours for a resonance mass 
of = 300 GeV. 


FIG. 1. Results for the observed exclusion likelihood, from the CMS ^ rr analysis [15], assuming 
a narrow resonance mass, m^, of 125 GeV (a) and 300 GeV (6). The solid (dashed) lines are obtained at 
^obs _ 2.28 (5.99) and indicate the approximate 68% (95%) C.L. allowed regions of a Higgs boson signal. 
The gray asterisk indicates the location of the global maximum of the likelihood. In (a) the yellow hollow 
diamond indicates the prediction of a Higgs boson at 125 GeV with SM signal strength. 


sum of the signal rates of the individual Higgs bosons in the cluster, 

(^{gg hj ^tt) = '^ a{gg hk) ■ BR(/ifc rr), (6) 

k 

a{gg bbhj rr) = ^ a{gg bbhk) ■ BR{hk rr). (7) 

k 

The cluster mass, mj, is determined by a signal strengths weighted mass average 

_ Efc [<^(gg ^fc) + crjgg bbhk)] • BR(/tfc rr) • nik 
Efc [^(gg hk) + a{gg bbhk)] • BR(/ifc rr) 

The sums in Eqs. run over all Higgs bosons hk combined in the cluster. In case 

there is no hj that fulhlls Eq. ([^ for a given hi, the cluster is formed solely by the Higgs 

boson hi- It should be noted that taking the incoherent sum of the contributions of the 
different Higgs bosons involves an approximation. While it is exact in the case of two 
different CV eigenstates, e.g. A and H in the MSSM, in general interference contributions 
can be important HSIES]. An extension of HiggsBounds that enables the implementation 
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of interference effects of nearby resonances in a generalized narrow-width approximation is 
currently under development, see also Ref. [20] , 

2. In the second step, the expected and observed likelihood values, 'i'lMdei’ 'respec¬ 

tively, for each Higgs cluster hj are evaluated from the experimental likelihood data grid. The 
likelihood is hrst evaluated for the rate values cr{gg —)• /17 —)• rr) and (7{gg —)• bbhj —)• rr), 
obtained through Eqs. in the mass-neighboring data slices, i.e. at the nearest grid 

mass values below and above mj, denoted by m_ and m_|_, respectively. The likelihood value 
at the predicted cluster mass mj is then obtained through linear interpolation: 

q(hi) = ■ i^+- + 9 +■ jrni - m.) ^ 

m+ — rri- 

Here q+/- denote the values of the test statistics obtained at the neighboring grid above or 
below the predicted mass m/ (we omitted the subscript ‘model’ for simplicity in Eq. Q). 
These are obtained, in each case, through bilinear interpolation within the two-dimensional 
likelihood planes of the provided CMS data. 

3. The steps 1. and 2. are repeated until all N neutral Higgs bosons have been evaluated as 
part of at least one Higgs cluster. 

4. Once all likelihoods have been evaluated, the most sensitive analysis application is deter¬ 
mined from the resulting expected likelihood, cluster is selected for which 

'i'nSdei(^^) niaximal. The observed exclusion likelihood, then used only for this 

cluster, and provides the hnal result. 

Hollowing this algorithm, the full likelihood from the CMS (j) ^ tt analysis for both the expected 
and observed exclusion can be directly obtained within HiggsBounds for any tested model. This 
is carried out via Fortran subroutines. For a technical documentation see Appendix |VI[ 

The use of this likelihood information is complementary to the other type of information con¬ 
tained in a full HiggsBounds application, which considers exclusion limits from many other Higgs 
searches from the LEP, Tevatron and LHC experiments. As an alternative to using the full like¬ 
lihood, we therefore also provide the option to reconstruct a limit at 95% C.L. and use this in 
the “standard” HiggsBounds operation. For clarity, we now repeat some elements of how this 
works [7]: In the statistical procedure, HiggsBounds hrst determines the most sensitive analysis 
to the model by picking the analysis application, for which the ratio between the model-predicted 
signal rate, ^predicted; over the expeeted upper limit on the signal rate, •S’^^Je^ted’ 

_ ‘S'predicted 

^expected = „95%CL ’ MUJ 

*^expected 

is maximized. After the most sensitive analysis has been determined, the model prediction is 
confronted with the observed exclusion limit of this particular analysis, The model is 
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considered to be excluded at the 95% C.L., if 


_ ‘S'pi'edicted 

’’observed = ^95%CL 


> 1 . 


( 11 ) 


observed 


For the CMS (/>—>■ rr analysis described above, and are a priori not known and 

need to be determined from the implemented likelihood distribution. In a numerical procedure, we 
therefore scale the model-predicted gg ^ cj) ^ tt and gg —)• bbcf) —)• rr rates with a universal factor 
/r until the obtained expected/observed likelihood values are equal to 5.99, corresponding to 

the two-dimensional 95% C.L. interval ^ The so-obtained scale factors, are then identified 

with the expected/observed 95% C.L. upper limits on the signal rate, respectively, which enter 
Eqs. (10) and (11). In this way, the likelihood-based results from the CMS ^ rr analysis can 


be incorporated in the standard HiggsBounds run. 


IV. VALIDATION 

Besides the nearly model-independent limits, CMS has also presented model-specific interpre¬ 
tations of their search results. This has been done for the MSSM, employing the benchmark 
scenarios proposed in Ref. [H] (see also Ref. [22] )• Here, we validate our likelihood implementation 
in HiggsBounds against the CMS results for three of these scenarios: The the light stop and 

the Iow-Mh scenarios (see Ref. [21] for details). The comparison of the reconstructed 95% C.L. ex¬ 
clusion line with the official CMS result provides a non-trivial test of our implementation: Firstly, 
it checks whether the exclusion likelihood agrees over a wide range of different compositions of 
the gluon fusion and b quark associated Higgs production rates obtained in the MSSM parameter 
space, which are mapped onto the two-dimensional likelihood grids (for fixed Higgs mass) in our 
reconstruction. Secondly, it tests whether our simple criterion of combining signal rates of Higgs 
bosons which have similar masses (overlapping within 20%) is a reasonable approximation. Thirdly, 
the validation also tests whether the results obtained from the statistical hypothesis test of a single 
narrow resonance model can be mapped reasonably well onto the full neutral Higgs spectrum of 
the MSSM (and beyond). 

Some deviations can be expected at the transition between regimes with different contribut¬ 
ing Higgs combinations. As explained above, the implementation in HiggsBounds is based on 
the CMS results for the single narrow resonance interpretation, and the contributions of different 
Higgs bosons of a considered model can only be combined if their mass differences are such that 
they would appear as a single resonance in the CMS search. In contrast, in the dedicated CMS 
analyses carried out in specific MSSM benchmark scenarios it was possible to properly combine 
the contributions from different Higgs bosons at any given mass constellation since these have been 
simulated and tested with their particular masses at every parameter point. Therefore, the dedi¬ 
cated CMS analysis is expected to have a higher sensitivity than the HiggsBounds implementation 


^ Technically, allowing for finite numerical precision we check for equality within < 1%. 
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(a) Distribution of the observed exclusion likelihood, (b) Map indicating the Higgs boson or cluster of Higgs 

'Zmssm! evaluated with HiggsBounds. The contours bosons with the highest sensitivity for a potential 

show the corresponding 95% C.L. exclusion limit exclusion that is used for the likelihood evaluation. 

{orange, solid) and the CMS result obtained from a 
dedicated analysis in this scenario [K] {green, dashed). 


FIG. 2. Exclusion likelihood evaluated with HiggsBounds in the (M^,tan/3) plane of the MSSM 
scenario. 


if multiple Higgs bosons with different masses each give a non-negligible contribution to the signal 
yield. Furthermore, due to the simple criterion used in HiggsBounds for including/excluding the 
contributions of additional Higgs bosons, the considered rates in HiggsBounds may change quite 
abruptly in a transition region, where the selection of the tested Higgs boson combination changes. 
The single resonance approximation is expected to work best when the signal can be described as a 
single resonance formed by one or several Higgs bosons, while contributions of other Higgs bosons 
besides those associated with the resonance are negligible. 

For predictions in the MSSM benchmark scenarios we employ the (M^,tan/3) grids of Higgs 
production cross sections and branching fractions for the MSSM benchmark scenarios provided 
by the LHC Higgs Cross Section Working Group (LHCHXSWG) [23] For the gg —> bh{h/H/A) 
production process we employ Santander-matching of the 4- and 5-flavor scheme (FS) cross sec¬ 
tions [52]. 


® The LHCHXSWG cross section and branching fraction grids for the MSSM benchmark scenarios are based on the 
following set of tools and calculations, that we list here for completeness: HIGLU [^, SusHi |25| . FeynHiggs | 26H31| . 
ggH@NNLO [32], HDECAY |3l|3l|, Prophecydf |35l[36j, bbh@NNLO(5FS) (37], bbh@NLO (4FS) [SEIES], ggH 
NLO massive [40], ggH NNLO for scalar Higgs [41] [4^, ggH NNLO for pseudoscalar Higgs [43] [44], EW corrections 
from light fermions [45l|46], (N)NLO (S)QCD corrections for hjHfA [471151] . 
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The results for the 


scenario in the plane are shown in Fig. ^ In Fig. 2(a) 


we show the distribution of the observed exclusion likelihood, (in color), as obtained from 

HiggsBounds. The corresponding 95% C.L. exclusion limit {orange, solid contour), which fulfills 
^MSSM ~ 5.99, is shown together with the CMS result obtained from a dedicated analysis in this 
benchmark scenario [T5| {green, dashed contour). As mentioned in Section [ITl the latter is based 
on the full combined 7 + 8 TeV dataset, whereas the exclusion information implemented in Higgs¬ 
Bounds is only based on the 8 TeV dataset. However, this fact is expected to lead to only minor 
differences in the excluded parameter regions. As can be seen, there is very good agreement between 
the exclusion limit reconstructed with HiggsBounds and the CMS result. Small deviations can be 
observed in the low Ma region, Ma ^ 150 GeV, where all three neutral MSSM Higgs bosons 
contribute substantially to the signal yield. Here, the result reconstructed with HiggsBounds 
excludes a slightly smaller area of parameter space. The HiggsBounds result can thus be regarded 
as a conservative estimate of the actual exclusion limit. 

In Fig. |2(b)| we display, for every parameter point in the (M^,tan/3) plane, the Higgs boson 
or combination of Higgs bosons (cluster) that has been selected to obtain the observed exclusion 
likelihood by the algorithm described in Section HI It can be seen that all three neutral Higgs 
bosons are combined in most of the parameter region with Ma ^ 160 — 170 GeV and tan/3 > 5, 
whereas at larger Ma values only the two heavier Higgs bosons, H and A, which are nearly mass 
degenerate, are combined to yield the most sensitive constraint. At low tan fd and large Ma values, 
however, the combined signal rate of the heavier Higgs bosons becomes so small that it is instead the 
light Higgs boson, with mass around 120 — 125 GeV, that is selected to give the observed exclusion 
likelihood. This is because its expected exclusion likelihood is larger than that obtained for H/A. 
The observed exclusion likelihood obtained for the light Higgs boson with mass around 125 GeV 
is non-zero because the best-fit point in the two-dimensional cross section grid at = 125 GeV 
is not identical with the SM prediction, cf. Fig. 1(a) This leads to the small, but non-zero 
values that are visible in Fig. 2(a) at large Ma and small tan/3. 

Next we look at the light stop benchmark scenario, for which the cross section predictions and 
their associated theoretical uncertainties have been discussed in detail in Ref. 1S3]- This scenario 
features a relatively low SUSY particle mass scale, MgusY = 500 GeV, and large stop mixing, 
Xt = 2 MgusY, leading to a lightest stop with a mass of ~ 325 GeV. This leads to a reduction 
of the gluon fusion cross section of the light Higgs by around 10 — 15% with respect to the SM 
prediction [21]. The results of applying our exclusion likelihood implementation in this scenario 
are shown in Fig. (with colors similar to Fig.[^. The agreement between the 95% G.L. exclusion 
contour obtained with HiggsBounds and the CMS result obtained from a dedicated analysis in 
this benchmark scenario [l5|, as displayed in Fig. 3(a), is very good for pseudoscalar Higgs masses 
Ma ^ 250 GeV. Similarly as in the scenario the reconstructed exclusion limit obtained from 

HiggsBounds is slightly weaker for lower Ma values than the CMS result from the analysis of the 
benchmark scenarios. As one can see in Fig. |3(b)[ the reconstructed likelihood in the Iow-Ma 

* Here, and in the following figures, we show as HiggsBounds result only the constraints obtained from the CMS 
(/ —>■ rr analysis and not from the full HiggsBounds application, where all currently implemented Higgs searches 
from LEP, the Tevatron and the LHC are taken into account. 
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Ma [GeV] 



100 200 300 400 500 600 

Ma [GeV] 


(a) Distribution of the observed exclusion likelihood, 
'Zmssm! evaluated with HiggsBounds. The contours 
show the corresponding 95% C.L. exclusion limit 
{orange, solid) and the CMS result obtained from a 
dedicated analysis in this scenario [K] {green, dashed). 


(b) Map indicating the Higgs boson or cluster of Higgs 
bosons with the highest sensitivity for a potential 
exclusion that is used for the likelihood evaluation. 


FIG. 3. Exclusion likelihood evaluated with HiggsBounds in the (ACi,tan/3) plane of the MSSM light stop 
scenario. 


region obtained from HiggsBounds is mainly based on a combination of the H and A signals, while 
in most part of this parameter region the light Higgs boson at 125 GeV is not covered by the 
20% mass overlap criterion used in HiggsBounds. In contrast, in the dedicated CMS analysis in 
this scenario also the contribution from the light Higgs boson h is properly combined with the 
contributions from the other neutral Higgs bosons. The latter analysis therefore has a slightly 
higher sensitivity in this region, which means that the exclusion bound that we find here is slightly 
conservative compared to the dedicated CMS result. In addition to the excluded parameter region 
at values of tan {3 > 5, the light stop scenario features an additional small excluded area at lower 
tan/? values, namely tan/? < 2, and Ma ~ 145 — 190 GeV. The exclusion contour evaluated with 
HiggsBounds matches very well with the CMS result also in this region, where gluon fusion is the 
dominant production mode. 

Finally, we test our implementation against the results obtained in the Iow-Mh scenario, where 
the heavy CP-even Higgs boson is interpreted as the discovered SM-like Higgs boson at around 
~ 125 GeV and the light CP-even Higgs is largely decoupled from the SM gauge bosons [M] , 
Unlike the other benchmark scenarios, which use Ma as a free parameter, this scenario is defined 
as a two-dimensional parameter plane in tan/? and the Higgsino mixing parameter pL. The mass of 
the pseudoscalar Higgs, Ma, is fixed to 110 GeV, which leads to a lightest Higgs mass that varies 
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low-m|_| scenario ‘^mssm 
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(a) Distribution of the observed exclusion likelihood, 
'Zmssm! evaluated with HiggsBounds. The contours 
show the corresponding 95% C.L. exclusion limit 
{orange, solid) and the CMS result obtained from a 
dedicated analysis in this scenario [K] {green, dashed) 
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(b) Map indicating the Higgs boson or cluster of Higgs 
bosons with the highest sensitivity for a potential 
exclusion that is used for the likelihood evaluation. 


FIG. 4. Exclusion likelihood evaluated with HiggsBounds in the (M^,tan/3) plane of the MSSM Iow-Mh 
scenario. 



mostly between ~ 80 GeV and ~ 105 GeV, but reaching even lower values at very low tan/3 and 
very high Since in the MSSM a low value of Ma implies also a light charged Higgs boson, this 
scenario served in particular as a benchmark for the LHC searches for charged Higgs bosons in top 
quark decays. In fact, the parameter space for this scenario in the MSSM is meanwhile essentially 
excluded [55l[57| (see also Ref. [58]). Nevertheless this benchmark scenario is still very useful for our 
validation since all three neutral Higgs bosons have similar masses and thus contribute non-trivially 
to the analysis. 


The comparison of the exclusion limits that have been reconstructed with HiggsBounds with the 
CMS results is shown in Fig.|^ It can be observed in Fig. 4(a) that there is rather good agreement 
between the exclusion limit obtained with HiggsBounds and the CMS result for ^ values up to 
/i ~ 2600 GeV. At this value of (depending on tan/3) the reconstructed exclusion limit develops 
an “edge”, and for higher /r values the reconstructed exclusion limit is significantly weaker than 
the one of the CMS result. The reason for this behavior is that at large /x the light Higgs mass 
becomes smaller than 88 GeV and is hence not combined with the heavier Higgs bosons A and H 
in HiggsBounds. This can be seen in Fig. |4(b)| In this parameter region, the tested signal rate 
is therefore signihcantly smaller than in the case of a full combination of h, H and A, and the 
resulting exclusion limit is accordingly weaker. In contrast, in the CMS analysis the signal yield 
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V Example application: “Alignment without decoupling' 


of the light Higgs h has been properly taken into account even for very low mass values, and the 
possibly decreasing signal efficiency is partially compensated by the increasing production cross 
section, leading to the significantly stronger exclusion at high ^ values obtained by CMSj^ 

In almost all the remaining parameter space, all three Higgs bosons are combined by Higgs- 
Bounds, as can be seen in Fig. |4(b)[ and the reconstructed exclusion contour resembles the CMS 
result for /i values below fj. ~ 2600 GeV. The slight deviations observed could result from mass- 
dependent selection efficiencies for the h and H signal yields, which cannot be accounted for in 
the HiggsBounds implementation since this information is not publicly available. Overall, even for 
this rather extreme scenario in the MSSM we find that the exclusion likelihood reconstructed with 
HiggsBounds approximates the results of a dedicated analysis reasonably well for large parts of 
the parameter space. 


V. EXAMPLE APPLICATION: “ALIGNMENT WITHOUT DECOUPLING” 


We now go beyond the validation with official CMS results and illustrate the usefulness of 
our exclusion likelihood implementation for another MSSM scenario. We consider here a scenario 
where the couplings of the light CP-even Higgs boson become SM-like for a certain range of tan (3 
values, independently of the masses of the remaining Higgs spectrum. The existence of this so- 
called alignment limit was first pointed out in Ref. |59) for the 2HDM. After the Higgs discovery 
this possibility has gained attention through a series of papers [HnHM] . see also the “r-phobic” 
benchmark scenario in Ref. [2T]. In the MSSM the alignment limit can be realized independently 
of the decoupling of the heavier Higgs states through a cancellation between tree-level and higher- 
order contributions in the Higgs sector. This cancellation can occur at relatively large values 
of tan /3 and jjL > Ms, with Ms = .ymbeing the stop mass scale. In the approximation 
tan/? ^ 1, and taking into account for simplicity only the dominant corrections at one loop, the 
alignment condition reads |64j 


tan j3 = 


"Uv|+45gj( 


2M| 


- 1 




6 Ml 


- 1 


( 12 ) 


Here, Mz and mt are the Z boson and top quark mass, respectively. M^ denotes the light CP-even 
Higgs boson mass in the above approximation, and v ~ 246 GeV. At is the trilinear soft-breaking 
term in the stop sector. 

Solutions of Eq. (12) with tan/3 > 0 exist if fj.At{A^ — 6M|) > 0. Typically, in order to achieve 
the correct Higgs mass Mh ~ 125 GeV for not too large values of the stop masses, the stop mixing 
is chosen in the region where the prediction for Mh is maximized, i.e. \Xt\ |At| ~ T-v/fiMs (at 
one-loop). Therefore, for fiAt > 0 {g-At < 0), the alignment condition has viable solutions for values 


® A better agreement in the large /j, parameter region could be obtained by increasing the mass overlap value of 20% 
in the criterion for forming Higgs boson combinations to a sufficiently high value. However, firstly, there is no 
strong physics motivation to choose a value well beyond the quoted mass resolution of ~ 20% of the experimental 
TT analysis. Secondly, values larger than 20% might lead to too aggressive exclusions in some scenarios. We 
therefore stick to the value of 20% as the default setting. 
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of \At\ that are slightly above (below) the value where the prediction for Mh is maximized. By 
increasing \^At/Mg\ it is possible to lower the tan/3 value at which alignment occurs. An MSSM 
benchmark scenario of this kind for BSM Higgs searches at the LHC has recently been proposed 
in Ref. [M] . 

Here, we investigate the benchmark scenario proposed in Ref. [Gjj; which is essentially a modifi¬ 
cation of the scenario [2T] to allow for alignment independent of decoupling. This so-called 

scenario is defined by the parameter values (in the on-shell scheme) 


Ml = 100 GeV, 


m, 


= mq = Mq, 


M 2 = 200 GeV, Ms = 1500 GeV, 

A^ = Aq = At, At/MQ = 2.45. (13) 

In contrast to the benchmark scenarios of Ref. |21] . the parameters and Mq are adjustable 
parameters in the scenario. For convenience, the slepton, sbottom and first and second 
generation squark soft-breaking mass parameters are set to Mq, however, these can easily be 
adjusted to higher values in order to avoid constraints from SUSY searches at the LHG as their 
influence on the Higgs phenomenology is negligible here. We follow the suggestion to set Mq to 
1 TeV per default and, if necessary, increase this value until a light Higgs mass of rrih > 123 GeV 
is obtained. In practice, this is only relevant at very low values of tan (5 in the benchmark scenario 
defined by Eq. (13). The parameter /x is then adjusted according to a chosen ratio ^/Mq. We focus 
here on the choice h/Mq = 3, implying rather large values of /r, where alignment independent of 
decoupling occurs at tan/3 ~ 10. 

The MSSM predictions are obtained using the public computer codes FeynHiggs-2.10.2 [26b 
[29ll31j for the Higgs masses, couplings and branching fractions, and SusHi-1.4.1 [25] for the gluon 
fusion and b quark associated production cross sections of the three neutral Higgs bosons. 

The numerical results for this benchmark scenario are displayed in Fig. The observed ex¬ 
clusion likelihood from the GMS (/>—)■ rr search as obtained from HiggsBounds is shown in color 
in Fig. 5(a), and the orange contour indicates the resulting observed 95% C.L. exclusion line. For 
comparison, the green contour shows the exclusion line obtained in Ref. [Mj using results from 


the same GMS analysis, however, following a more simplistic approach]^ As can be seen from the 
figure, the more advanced implementation of the observed exclusion likelihood in HiggsBounds 
leads to a somewhat stronger 95% G.L. exclusion limit over most of the parameter space. The 
relative behavior seems to be different in the ti threshold region, Ma ~ 2mt ~ 345 GeV, where in 
particular the gg ^ A cross section is enhanced. However, the approximations made in Ref. [Mj 
appear to be least reliable in this region. The gray dotted line shows the exclusion limit obtained 
by GMS for the scenario. As discussed in Ref. |2T], the excluded regions in the benchmark 

scenarios are significantly affected if decay modes of the heavy Higgs bosons H and A into SUSY 
particles are kinematically open and unsuppressed. The presence of such decay modes leads to 


® The limit in Ref. has been obtained by “reverse-engineering” an inclusive [<y{gg —^ </>)—>■ hhcf)] x BR(</) -A 
tt) limit from the CMS results for the scenario with g = 200 GeV m, and applying this cross section 

limit to the given alignment benchmark scenario. In particular, this approximation does not take into account 
the sensitivity of the limit on the individual contributions from the gluon fusion and b quark associated Higgs 
production processes. 
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Ma [GeV] 


(a) Likelihood from hjHjA —>■ rr exclusion. 


(b) Likelihood from Higgs signal rates. 


FIG. 5. Constraints from LHC Higgs searches in the alignment benchmark scenario (with /r = 3Mq): 
(a) Distribution of the exclusion likelihood from the CMS (j) ^ tt search and observed 95% C.L. exclusion 
line as obtained from HiggsBounds. For comparison, also the corresponding 95% C.L. exclusion line given in 
Ref. |Mj {green, solid) and the 95% C.L. exclusion line in the scenario with /r = 200 GeV obtained 

from HiggsBounds {gray, dashed) are shown. (&) Likelihood distribution, A^^g, obtained from testing the 
signal rates of the light Higgs boson h against a combination of Higgs rate measurements from the Tevatron 
and LHC experiments, obtained with HiggsSignals. The minimal is found at the gray asterisk. 


a sizable reduction of the H/A -A tt branching fractions and therefore to a smaller excluded re¬ 
gion. In the alignment scenario [i is very large, leading to a negligible Higgsino component in the 
light neutralinos and chargino. The branching fractions for the Higgs decays to neutralinos and 
charginos are therefore essentially absent. In addition, the heavy Higgs decays to gauge bosons, 
H -A W~^W~ and H -A ZZ, are also suppressed, as the responsible coupling oc cos(/3 — a) vanishes 
in the alignment limit. As a result, the H/A -A tt branching fraction is significantly higher in the 
alignment scenario than in the 771 ™°'^"’" scenario, which leads to a much larger excluded region in 
the alignment scenario, see also the discussion in Ref. |64j . 


In order to illustrate the complementarity between the constraints from the CMS (/>—>■ rr search 
and the constraints obtained from the signal rate measurements of the discovered Higgs boson, we 
show in Fig. 5(b) the likelihood distribution, Axhsj obtained from a test of the light Higgs boson 
signal rates against a combination of the latest rate measurements from the LHC |65H73] and the 


Tevatron [HITS], using the public computer code HiggsSignals-1.3.0 [8] (see also Refs. [UlIZS])- 
The 95% C.L. preferred region lies within the orange contours in Fig. 5(b) It is given by the x^ 
difference with respect to the minimal x^ value (located in the alignment region and indicated as 
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FIG. 6. Combination of constraints from the CMS (/)—>■ tt search and the latest Higgs rate measurements 
in the MSSM alignment scenario (with /i = 3Mq): The global function, AxtotJ based on the likelihoods 
provided by HiggsBounds and HiggsSignals, is shown in color; The contours indicate the Icr, 2cr and 3(t 
allowed regions. 


gray asterisk in Fig. 5(b)[ ), Axhs = Ahs 
becomes independent of Ma at around tan j3 


Xhs,: 


^in < 5.99. It can be seen that the distribution 
10, indicating that the couplings of the light Higgs 


become SM-like independently of the decoupling of the heavier Higgs states. 


Since we now have the exclusion likelihood from the CMS (/> —)> rr search available, we 

can perform a statistical combination with the constraints from the Higgs rate measurements by 
constructing the global function = (Z^ggjyt + X^S ■ resulting Axtot distributiorj^ is shown 
in Fig. The constraints from the rr searches at the LHC are highly complementary to 
the rate measurements, since they are particularly sensitive at higher values of tan /3 where the 
production process gg —?■ bbcj) is enhanced. In the scenario with g = 3Mq, the combination 
of both constraints yields a lower limit of Ma ^ 350 GeV at the 95% C.L. Thus, alignment of 
the light Higgs boson occurring without the simultaneous decoupling of the heavier Higgs states is 
ruled out for this scenario. The alignment without decoupling limit can be pushed to lower values of 
tan/3 in this scenario, where the constraints from the (/i —)■ rr searches are less significant, only by 
choosing even more extreme values of gAt/Mg, which potentially leads to problems with vacuum 
stability ITT]. 


^ Again, Axtot is the x^ difference with respect to the minimal x^ value (obtained at Ma = 500 GeV, tan/3 = 4, 
i.e. in the lower right corner of Fig. j^, now based on the global likelihood Xtot- 
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VI. CONCLUSIONS 

LHC searches for non-standard Higgs bosons decaying into tan lepton pairs constitute a sensitive 
experimental probe for BSM physics. Recently, the CMS collaboration published the likelihood 
information for their Higgs boson searches in the final state m- The likelihood is given 

as a function of the two relevant Higgs production channels, gluon fusion and b quark associated 
production, for various mass values of the narrow resonance assumed for the signal model. In this 
paper we have shown how this experimental information can be utilized to test large classes of 
theoretical models. In particular, we have developed a simple algorithm that maps an arbitrary 
model with multiple neutral Higgs bosons onto a model with a single narrow resonance, for which 
the corresponding exclusion likelihood from the CMS search can be determined. We have described 
the inclusion of this method into the new version of the publicly available Fortran code HiggsBounds 
(version 4.2.0 and higher). For nearly any model under consideration, HiggsBounds provides an 
evaluation of the exclusion likelihood for a model parameter point based on the information from 
Ref. [15| . Similarly, if requested, HiggsBounds can also perform a test of whether or not a given 
parameter point is excluded at the 95% C.L. based on all available searches, including the new 
T~^T~ result. The approach to test BSM models with exclusion limits is complementary to testing 
the compatibility of a given model with the observed Higgs signal (and possible future signals 
of additional Higgs bosons). The latter kind of information is contained in the sister program 
HiggsSignals, and both programs can be used together in order to obtain the combined likelihood 
information from the search limits and the Higgs rate measurements. Both codes are available at 
http://higgsbounds.hepforge.org, 

We have validated onr implementation of the search resnlts into HiggsBounds by compar¬ 

ing the 95% C.L. exclusion contours obtained with HiggsBounds with the ones obtained by CMS 
from dedicated analyses in three Higgs benchmark scenarios |21) in the MSSM. We found very good 
agreement in the parameter regions where the sensitivity of the search is dominated by a single 
combination of Higgs bosons that can be identified with a single narrow resonance assuming an 
experimental mass resolution of 20%. As expected, the largest but still relatively small deviations 
occur in parameter regions where all neutral MSSM Higgs bosons are relatively close in mass and 
contribute comparably to the signal yield. 

As an application, we have discussed the combined constraints of the rr search and the rate 
measurements of the SM-like Higgs at 125 GeV in a recently proposed MSSM benchmark scenario, 
where the lightest Higgs boson obtains SM-like couplings independently of the decoupling of the 
heavier Higgs states. Here we combined the analysis of the rate measurements for the Higgs 
signal, evaluated with HiggsSignals, with the exclusion likelihood from the non-observation in the 
T~^T~ search channel, evaluated with HiggsBounds. We have shown that the combined information 
yields very significant constraints on the available parameter space in this scenario and in fact 
disfavors the “alignment without decoupling region” in the studied benchmark model. 

We encourage ATLAS and CMS to continue providing their search results including the relevant 
likelihood information. This will greatly facilitate the application of the search results for testing 
BSM models. 
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USER GUIDE: HOW TO OBTAIN THE EXGLUSION LIKELIHOOD WITH HIGGSBOUNDS 


The exclusion likelihood information for the GMS (?!)—>■ rr analysis |15j is implemented in 
HiggsBounds from version 4.2 on. As described in Section [ITI| this information is used in a standard 
HiggsBounds run to reconstruct the expected and observed 95% C.L. exclusion limit, which is then 
considered alongside all other available Higgs search limits in the full HiggsBounds application. 
This leads to the global information whether the tested parameter point is allowed or excluded 
at the 95% C.L. Beyond this information the value for the exclusion likelihood for the model 
parameter point under investigation, gmodeU can also be obtained directly via HiggsBounds Fortran 
subroutines, enabling the user to incorporate this information e.g. in a global parameter fit. In the 
following we document the relevant subroutines that make this information accessible. 


The main routine that runs the algorithm presented in Section HI to obtain the exclusion 
likelihood is: 


/-N 

HiggsBounds_get_likelihood( int analysisID, int Hindex, int nc, int cbin, dhle M, 

dhle llh, char(*) obspred) 

The (mandatory) input argument|^ analysisID specihes the analysis for which the likelihood 
should be obtained. At the moment, the CMS 0 —>■ rr analysis based on the full 8 TeV dataset 
(analysisID = 3316) is the only available likelihood, but the framework is easily extendable for 
future experimental results. The output values provide information about the selected Higgs boson 
combination (or Higgs cluster): 


Hindex gives the index i of the Higgs boson hi, which provided the initial seed to form the 
dominant Higgs cluster (cf. Section HI item 1), 


- nc gives the number of Higgs bosons contained in the combination. 


* Here, and in the following, input arguments and optional input arguments are highlighted in dark blue and green, 
respectively. The remaining arguments are output values. 
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- cbin is a binary code (bitmask) for the identifiers of the participating Higgs bosons. The 
binary code is given by summing over for all involved Higgs bosonsj^ For example, in 

the MSSM (with common indexing hi = /i, /12 = H, = A), the combination H+A would 
give cbin = 6, whereas a cluster formed only by the light Higgs h gives cbin = 1. 

The output value M gives the averaged mass value, calculated according to Eq. Q, at which the 
likelihood value has been evaluated. The computed value of the likelihood, gmodeb is returned 
as llh. The final argument is an optional input, obspred, which takes a string value that can 
be either ‘obs’ or ‘pred’, specifying whether the observed or expected (or predicted) exclusion 
likelihood should be evaluated, respectively. The default behavior if this argument is not provided 
is that the routine returns the observed likelihood, following the algorithm described in Section [111} 
In addition to the main subroutine we provide the following two auxiliary routines, which may 
be helpful to understand the obtained results. 


HiggsBounds_get_likelihood_for_comb( int analysisID, int cbin.in, int Hindex, 

int nc, int cbin, dble M, dble llh, 
char(*) obspred) 

This routine evaluates the likelihood for a specific selection of Higgs bosons that should be consid- 
ered for the test. Higgs bosons that are not available for a possible formation of a Higgs cluster are 
specified with the input parameter cbin.in, which is a binary code following the same convention 
as cbin above. The remaining arguments are the same as above, with the only exception that 
obspred is now a mandatory input parameter. Among the available Higgs bosons, the routine 
selects the Higgs combination with the maximal likelihood value and provides the corresponding 
results. 


HiggsBounds_get_likelihood_for_Higgs( int analysisID, int cbin_in, int Hindex, 

int nc, int cbin, dble M, dble llh, 
char(*) obspred) 

I’his auxilliary subroutine works in a similar way as above. However, the additional input argument 
Hindex forces the routine to consider only Higgs clusters that contain the specified Higgs boson hi. 

In global parameter fits, where both the conventional HiggsBounds output (95% C.L. exclusion) 
as well as the likelihood information is used, it is often convenient to deactivate specific analyses 
during the standard HiggsBounds run, since these are better described by the likelihood informa¬ 
tion. In particular, the 95% C.L. limits from previous BSM </> —?• rr searches can be deactivated 
if instead the CMS ^ rr exclusion likelihood is used. In order to do so, we provide two new 
subroutines, contained in the Fortran module ‘ channels L 


The indexing of the Higgs bosons is identical to the ordering in which the user chooses to specify the theoretical 
input for HiggsBounds [lEllT]. 
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cUialysisID 

Experiment Luminosity and GM-Energy 

Additional notes 

Reference 

3316 

GMS 

19.7 fb”" at 8 TeV 

using — 21nL reconstruction 

m 

2014049 

ATLAS 

19.5 - 20.3 fb"^ at 8 TeV 

profiled limit on gg — >■ bbf process 

US] 

20140492 

ATLAS 

19.5 - 20.3 fb"^ at 8 TeV 

profiled limit on gg ^ f process 

m 


TABLE I. Implemented 95% C.L. exclusion limits from LHC searches for BSM Higgs bosons with rr final 
states in HiggsBounds-4.2. The analysisID is used as a unique identifier for an individual analysis and 
can be used to deactivate/activate them in HiggsBounds (see text). 


— 

HiggsBounds_deactivate_aiialyses ( int (:) aiialysisID_list) 

«- 

This routine should be called before the subroutine run_HiggsBounds in order to deactivate the 
analyses specified by the integer array analysisID_list. For convenience, the analysis identifiers 
of the currently implemented LHC (j) ^ tt searches are given in Tab. |I| 


HiggsBounds_activate_all_analyses 0 

This subroutine can be used at any time to re-activate all previously deactivated analyses for the 
succeeding HiggsBounds run. 

In order to demonstrate the use of these subroutines, we provide the example program 
HBwithLHClikelihood, included in the /example.programs/ directory of the HiggsBounds distri¬ 
bution. This program shows how to obtain the observed exclusion likelihood from the CMS (/> —)■ rr 
analysis in the MSSM scenario, such that the user should be able to directly reproduce Fig.j^ 
The example can be compiled by calling ‘make HBwithLHClikelihood’ in the HiggsBounds main 
folder, and run from the example.programs folder by calling ‘ ./HBwithLHClikelihood’. Follow¬ 
ing a successful run, the gnuplot scripts ‘plot_mhmax_llh.gnu’ and ‘plot_mhmax_llh_comb.gnu’ 
in the same folder then reproduce Fig. 
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